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ABSTRACT 

We study the evolution of the bias factor b and the mass-galaxy correlation coefficient r in a simple 
analytic model for galaxy formation and the gravitational growth of clustering. The model shows that b 
and r can be strongly time-dependent, but tend to approach unity even if galaxy formation never ends as 
the gravitational growth of clustering debiases the older galaxies. The presence of random fluctuations 
in the sites of galaxy formation relative to the mass distribution can cause large and rapidly falling bias 
values at high redshift. 

Subject headings: galaxies: statistics — large-scale structure of universe 



1. INTRODUCTION 

The relative distribution of galaxies and mass is of in- 
creasing concern in cosmology. Constraints on cosmolog- 
ical parameters from galaxy surveys are only as accurate 
as our understanding of bias. Furthermore, if the den- 
sity parameter in mass that is capable of clustering were 
shown to be O = 0.25 ± 0.15 then most dynamical anal- 
yses of the motions of galaxies would be consistent with 
the assumption that galaxies trace mass, and the challenge 
would be to explain why simulations of the adiabatic cold 
dark matter (CDM) model for structure formation typ- 
ically indicate galaxies are more strongly clustered than 
mass (Jenkins et al. 1998). If £1 = 1, galaxies do not trace 
mass and it would be puzzling that the mass autocorrela- 
tion function in CDM simulations is a much poorer approx- 
imation to a power law on scales 10 kpc < hr < 10 Mpc 
than is the galaxy two-point correlation function, and that 
the mass function shows a much more pronounced evolu- 
tion of shape back to z ~ 1 (Jenkins et al. 1998). Major 
surveys in progress of galaxies and the cosmic microwave 
background (the CMB) will advance our understanding of 
such issues. It may be useful to supplement these observa- 
tions with analytic illustrations of how biasing can evolve 
as mass clustering grows and gravity draws together the 
galaxies with the mass. The purpose of this Letter is to 
present a simple analytic model for the dynamical evolu- 
tion of the relative distribution of galaxies and mass. 

We assume that at formation a galaxy may be assigned a 
near permanent observational tag. The far infrared lumi- 
nosity would not do, for IRAS galaxies avoid dense regions, 
an example of biasing due to environment (e.g., Strauss & 
Willick 1995). The spheroid luminosity may be an ade- 
quate tag, for the spheroid star populations in more lumi- 
nous galaxies are thought to be old and slowly evolving. 
We let pz(r) be the density of this luminous matter and 
p(r) be the total mass density, both smoothed on some ap- 
propriate scale R. The corresponding density fluctuations 
are S = p/(p) — l and 5i = pi /(pi) — 1. In a commonly used 
model Si = b8, where the constant b is the bias factor. We 
adopt the more general and possibly more realistic statis- 



tical representation of Dekel & Lahav (1998; Dekel 1997 
§5.5; Lahav 1996 §3.1) where 8 and 5i are treated as sta- 
tionary (translationally invariant) random processes that 
may be grouped in the two-dimensional vector 
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The mean (x) vanishes by definition. We write the covari- 
ance matrix as 
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where a = C 1 { is the rms mass fluctuation on the scale 
R, b = (C22/C11) 1 / 2 is the bias factor (the ratio of lumi- 
nous and total fluctuations), and r = C12/ '(C11C22) 1 / 2 is 
the dimensionless correlation coefficient between the dis- 
tributions of mass and galaxies. If the smoothing scale 
R is large enough that x can be modeled as a bivariate 
Gaussian random variable, then C contains all the statis- 
tical information about x. In principle this representation 
is not affected by galaxy merging, if the process does not 
add much to the star population that tags galaxies. In 
practice many of the surveys in progress will use galaxy 
counts, but it is thought that there has not been sub- 
stantial merging of L > i» galaxies since redshift z = 1. 
Pen (1998) has shown that the three quantities a, b and 
r can be measured from redshift space distortions (Kaiser 
1987; reviewed by Hamilton 1997). High redshift surveys 
in progress (see e.g. Giavalisco et al. 1998; Yee et al. 1998 
and the review in Moscardini et al. 1998) open the possi- 
bility of measuring the time evolution of b and r. Thus it 
seems timely to explore physical models for the evolution 
of these second moment measures of biasing. 

For studying the scale-dependence of bias, it is useful to 
work with x, the Fourier transform of x, which satisfies 
(x(k)x(k') f ) = (27r) 3 <5 D (k - k')C(k), 
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Here P is the power spectrum of the mass distribution, Pi 
is the power spectrum of the light that traces the galaxies, 
and P x is the cross spectrum. The last expression defines 
the analogs of the bias functions in equation (|J). We have 
entered vector arguments to take account of application 
in redshift space in fields of small angular width. The fol- 
lowing analysis applies equally well to b and r as to b' and 
r'. 

Most previous analyses of biasing models have focused 
on the static and local or non-local relation between Si 
and 5 as galaxies form. Recent reviews are given by Mann 
et al. (1998), Croft etal (1998) and Scherrer & Weinberg 
(1998); the latter appears to be the first to make static 
model predictions for r as well as b. But these measures 
are functions of time. Thus even if galaxies initially were 
uncorrelated with the mass (r — 0), they would gradually 
become correlated as gravity draws them towards over- 
dense regions, and one might expect this process to drive 
b and r toward unity. Fry (1996) has demonstrated this 
explicitly for the special case r = 1, and similar conclu- 
sions have been found in numerical simulations of dark 
matter halo clustering (e.g., Mo & White 1996; Matarrese 
et al. 1997; Bagla 1998; Catelan 1998ab; Porciani 1998; 
Wechsler et al. 1998). 

We limit our discussion to the linear perturbation theory 
regime \6\ <C 1, that is, to large scales R. In §2 we derive 
a general expression for the time-evolution of C that ap- 
plies once galaxy formation has stopped. We generalize to 
ongoing galaxy formation in §3, and present conclusions 
in §4. 

2. AFTER THE GALAXY FORMATION EPOCH 

We assume galaxy formation converts part of the mass 
into a near permanent luminous form without affecting its 
position or velocity on large scales, and that this lumi- 
nous form behaves as test particles that are fair tracers of 
the large-scale velocity field (but not necessarily the mass 
density) . Since galaxies and mass are assumed to have the 
same bulk peculiar velocity v(r), the contrasts Si and 5 sat- 
isfy the same linear continuity equation, &i « 8 w —V • v, 
in the absence of galaxy formation. This gives 
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The subscript indicates the present value, and D is the 
growth factor of the mass density contrast in linear per- 
turbation theory (see e.g. Peebles 1980) normalized to the 
present value D(oq) = 1. The covariance matrix therefore 
evolves to 

C = MCM*. (5) 
Equations (||) to (||) give 



cr = cr/D, 

b a = [(1 - Df + 2D{1 - D)br 
r = [(1 -D) + Dbr]/b , 
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The first equation simply states that the total mass fluc- 
tuations have grown as D(t). The second two equations 
show that the situation tends to grow simpler, b and r ap- 
proaching unity regardless of their initial values. This is 
illustrated in Figure 1. Inverting equations ^ to (0) gives 
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2(l-D)b r + b 2 ] 1 / 2 /D, 
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This tells us what a, b and r must have been in the past 
to produce the present values. In an Einstein-de Sit- 
ter universe (Q m = 1 and Q\ — 0, as in the standard 
SCDM cold dark matter model), the linear growth factor 
is D = a/ao = (1 + z)^ 1 and these equations are 

(12) 
(13) 
(14) 



z) V , 



a=(l 

b=[z 2 - 2z(l + z)b Q r + (1 + zfb 2 ] 1 / 2 , 
r=[(l + z)b Q r -z]/b. 

For reference, these SCDM redshifts z are given at the top 
of Figure 1 . 

The value r = 1 is a fixed point: if mass and galaxies 
are perfectly correlated at one time this remains true for 
all time. For r = 1, equation (|?|) is 

(b -l) = (b- 1)D, (15) 

which for SCDM gives b = bo + (bo — l)z, the special case 
derived by Fry (1996) and Mo et al. (1997). This is the top 
curve in each set in Figure 1A. One also sees that if r <~ 0, 
b tends to decrease even if there is no initial bias. The 
correlation r grows monotonically with time in all cases. 
The larger b, the slower the approach of r to unity. 
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FIG. 1 — The evolution of bias (A) and correlation (B) is shown 
for 15 models that at Einstein-de Sitter redshift z = 5 have bias 
b = 0.5, 1, and 2 and correlation coefficients r = 0, 0.25, 0.5, 0.75, 
and 1. In the top panel, r increases upward in each quintuplet of 
lines. In the bottom panel, bo increases downward in each triplet of 
lines. The present epoch is at D~ 1 = 1, and the evolution is extrap- 
olated to a large future increase in the density contrast, D^ 1 — » 0. 
The redshift scale at the top of the figure assumes the Einstein-de Sit- 
ter model. 
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3. ONGOING GALAXY FORMATION 

We consider now the more general case where galax- 
ies form over a substantial range of redshift. Here there 
is a competition between gravitational instability, which 
pushes b and r toward unity, and the formation of new 
galaxies, which can be biased or poorly correlated with 
the mass distribution. We model the time evolution of 
the density of the galaxies, as measured by the luminosity 
tracer pi, as 



pi + Pl V ■ v = g(t) [1 + M + S±] 



(16) 



in comoving coordinates. The galaxy formation rate per 
unit volume is g(t). The dimcnsionless factor in brackets 
models the sum of galaxy formation determined by the lo- 
cal mass density and a random component S± (r, t) . The 
latter has zero mean ((8±) = 0) and is uncorrelated with 
the mass density ((88 ±) = 0). The deterministic part is 
represented as the first two terms in a Taylor series expan- 
sion of the galaxy formation rate as a function of the local 
mass density, as in Coles (1993) and Fry & Gaztanaga 
(1993). The time-dependent parameter &» in this expan- 
sion is the "bias at birth" . 

As in the previous section, the streaming velocity of 
galaxies and mass is related to the time-evolution of the 
mass contrast, — V • v = 8 = D8q. The space average of 
equation ( |i~6| ) is to leading order (pi) = g(t), so the mean 
galaxy density is 



(17) 



(fn{t)) = G(t)= / g(t')dt'. 
Jo 

Thus we can rewrite equation ( |l6| ) as 

pi = g + (KDg + DG)8 + gS X) 



(18) 



and we integrate this to obtain pi . For any time-dependent 
quantity f(t) we define the time average weighted by the 
galaxy formation rate g by 
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With this notation and the relation DG = (DG)' — Dg, 
equation (|18|) gives the galaxy contrast 
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Pi 
(pi) 

= i ^ 



1 = cS+(8 ± ) t , 



((K - l)D) t 



D 



Then the covariance matrix is 
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in the notation of the last term in equation ( |l9| ) . Combin- 
ing equation (0) with equation (]22|) yields 



(24) 



[c 2 + S 2 /a 2 ] 1 ' 2 , r = c/b. 



Let us first consider some simple cases. If there is no 
random contribution to galaxy formation, 8j_ = 0, then 
s 2 = 0, r = 1, and b = c. If r = 1 and the bias at 
birth, is time- independent, then equation ( pl| ) reduces 
to (bo — 1) = (6* — 1)(-D)t- Here the deviation of the bias 
from unity is suppressed by the average growth factor be- 
tween the galaxy formation epoch and today. This is just 
equation (fLq) averaged over the galaxy formation history. 
In the the Einstein-de Sitter model this bias suppression 
factor (D) t is ((1 + -z)" 1 ^. 



Table 1 - Resulting bo and ro for various models, for a constant 
galaxy formation rate g(z) between redshifts z on and z //- 



Model 


n 


A 




Z off 


b, 


s* 


N 




Ml 


i 





5 


5 


2 





1 


1.17 


1 


M2 


0.3 


0.7 


5 


5 


2 





1 


1.21 


1 


M3 


0.3 





5 


5 


2 





1 


1.30 


1 


M4 


1 





5 





2 





1 


1.36 


1 


M5 


1 





1 





2 





1 


1.69 


1 


M6 


1 





3 


1 


2 





1 


1.35 


1 


M7 


1 





5 


2 


2 





1 


1.23 


1 


M8 


1 





7 


3 


2 





1 


1.17 


1 


M9 


1 





5 


2 


0.5 





1 


0.88 


1 


M10 


1 





5 


2 


1 





1 


1 


1 


Mil 


1 





5 


2 


1 


0.3 


1 


1.08 


.93 


M12 


1 





5 


2 


1 


0.3 


10 


1.01 


.99 


M13 


1 





5 


2 


2 


0.3 


1 


1.29 


.95 


M14 


1 





5 


2 


2 


0.3 


10 


1.24 


.99 



Redshift z 
3 2 1 




4 3 2 

Inverse growth factor 1 D 1 

FIG. 2 — The time-evolution of bias (A) and correlation (B) for 
some of the models in Table 1. 



The randomness contribution s 2 to the correlation r de- 
pends on (S±(t')8±(t")) , which is a measure of the memory 
of fluctuations 8± at the same comoving position at dif- 
ferent times. If memory were complete, 8j_ a function of 
position alone, then s = (i5^)^ 2 . In the limit of a short co- 
herence time, (8±(G')8_ L (G")) oc 8(G'-G"), equation tfxfy 
gives s oc G -1 / 2 , just as the average of N independent 
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fluctuations is proportional to N 1 ^ 2 . The simple two- 
parameter model 



S*0"Q 



1 + (N 



(25) 



incorporates both of these extremes as special cases (N = 
1 and N = oo, respectively), and since this aspect of the 
galaxy formation process is still so poorly understood a 
more complicated model for s does not yet seem war- 
ranted. The parameter TV can be interpreted as the ef- 
fective number of uncorrelated galaxy formation epochs. 
The current rms mass contrast is <7o, and the other model 
parameter is the relative normalization s*. 

Examples of these relations are listed in Table 1 and 
plotted in Figure 2. The growth factors D for the low 
density models are computed as in Carroll et al. (1992). 
These models show that the biasing parameters r and b 
can vary quite rapidly near the start of galaxy formation, 
and that r and b tend to unity even if galaxy formation 
never ends, as in models 4 and 5, because the growth of 
gravitational clustering debiases the growing number of 
older galaxies. After galaxy formation terminates b and r 
approach unity more rapidly, evolving as in Figure 1. A 
random component S± produces a large early value of b, 
but b and r approach unity quite rapidly if N 3> 1, as the 
fluctuations from large numbers of random events average 
down to insignificant levels. 

4. CONCLUSIONS 

We have computed the evolution of the bias factor 6 and 
the mass-galaxy correlation coefficient r in a simple ana- 
lytic model represented by a galaxy formation history g(z), 
the parameters in a representation of biased galaxy forma- 
tion (6*, a*, and N), and the parameters for the cosmolog- 
ical model (fi and A). Galaxy formation could be affected 
by explosions, as from quasars, at positions unrelated to 
the local mass density (e.g., Dekel & Rees 1987; Babul & 
White 1991). This effect is represented by the term Sj_. 



The coherence time of 5±, as measured by N, would de- 
pend on how long individual quasars last and how long 
the intergalactic medium "remembers" their effects. The 
complications of luminosity evolution and galaxy mergers 
are incorporated into our formalism by defining g(z) to be 
the formation history of the matter that is luminous at the 
rcdshift z' at which we wish to evaluate b and r. 

If galaxies were assembled at z < 1 there could be a 
large and observationally significant evolution of the bi- 
asing functions b and r at low redshift, as illustrated by 
the steep initial slopes in Figure 2. On the other hand, if 
galaxy positions were assigned at relatively high redshift 
the effect of biasing at birth would be strongly suppressed. 
Decreasing the mass density fi decreases this effect, be- 
cause the growth of density perturbations at low redshifts 
is slower, but this likely is at least partially canceled by 
shift of the formation epoch to larger redshifts. In fact, if 
g, 6*, and S± were functions of D alone, debiasing would 
be independent of the cosmology. 

The new galaxy redshift surveys, including the two de- 
gree field (2dF) survey and the Sloan Digital Sky Sur- 
vey, have the potential to measure key cosmological pa- 
rameters with great accuracy, both alone (Tegmark 1997; 
Goldberg & Strauss 1998), and when combined with cos- 
mic microwave background experiments (Hu et al. 1997). 
To achieve this, however, biasing and its evolution must 
be understood to a comparable accuracy. Our formalism 
deals with the second moments of the mass and galaxy dis- 
tributions in a universe that is statistically homogeneous 
(eq. H ) . The trend to debiasing of the functions b and r in 
this representation appears to be a generic feature of the 
gravitational instability scenario for structure formation, 
a prediction that may be observationally tested in the near 
future as galaxy redshift surveys improve. 

We thank Ben Bromley and Ue-Li Pen for useful dis- 
cussions. Support for this work was provided by NASA 
though grant NAG5-6034 and a Hubble Fellowship, HF- 
01084.01-96A, awarded by STScI, which is operated by 
AURA, Inc. under NASA contract NAS5-26555, and at 
Princeton University by the National Science Foundation. 
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